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Q . Abstract. The empirically determined universal power-law shape of X-ray 

[ luminosity function of high mass X-ray binaries in galaxies is explained by 
fundamental mass-luminosity and mass-radius relations for massive stars. 
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lO ' The unique possibilities offerred by modern X-ray space observatories Chan- 

^ ! dra and XMM-Newton make it possible to observe point-like X-ray sources 

fS| ! in nearby galaxies which allows studies of formation and evolution of the 

O I population of close binary systems as a whole. Recently, Grimm, Gilfanov 

1^ • and Sunyaev (2002, 2003) reported that the differential luminosity function 
of point-like X-ray sources in nearby galaxies has a universal power-law shape 
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^ : ^ cx SFR X L;" 

a ■ dLx 

in the wide range of X-ray luminosities ~ 10^^ — 10^" erg/ s, where SFR means 
r\ [ the star formation rate in a specific galaxy, a = 1.61 ± 0.12. 
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In this Letter we show that the empirically found power-law form of the 
X-ray luminosity function naturally appears during accretion of matter in 
close binary systems. In high-mass X-ray binaries, the compact star accretes 
matter from stellar wind from an early-type optical companion. Observations 
suggest (De Jager, 1980) that stellar wind mass loss rate from massive (> 
several Mq) main sequence stars can be parametrized as 

oc — (1) 

where L is the stellar luminosity, Voo is the stellar wind velocity at the infinity. 
With a good accuracy the latter is proportional to the escape velocity near 
stellar photosphere foo ~ Svp. Formula (|l|) simply refelcts the radiation 
pressure driven wind. 
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It is well known that for stationary stars on the main sequence power- 
law mass-luminosity (M-L) and mass-radius (M-R) relations hold (see e.g. 
Zel'dovich, Blinnikov, Shakura 1981). The M-L relation follows from the 
diffusive character of the heat transfer from stellar interiors 

L oc l/ztpVT^ (2) 

(k is the opacity, p is the density, T is the temperature) with account of 
the virial relation T cx M/R for stationary stars. For sufficiently massive 
stars Tompson scattering dominates in the opacity and L oc (Eddington, 
1926). The empirically determined power-law index for early-type stars is 
close to the theoretically expected value 3, with slightly decreasing down to 
2.7 for very massive stars (De Jager, 1980). 

The M-R relation can be obtained by using the energy release per gramm in 
the stellar core due to thermonuclear reactions 

e ~ pT^' (3) 

where Ze = (9 log < avo > /d\ogT is the Zel'dovich number showing the 
temperature dependence of thermonuclear reactions (see Blinnikov 2000) and 
radiation diffusion equation (^): i? oc M^*^"^/^*^"*"^. For example, for the pp- 
cycle Ze ~ 3. ..10, for the CNO-cycle Ze ~ 10. ..27, depending on the central 
temperature. The empirically found mass-radius relation for main sequence 
stars is i? oc M^^^ with some decrease in the power-law index for very massive 
stars. 

The accretion rate onto a compact star from stellar wind is determined 
by its orbital velocity (fort), the masses by the binary system companions 
(Mj; and Mq for the compact and optical star, respectively) and by the stellar 
wind velocity (Lipunov 1992). In the case of fast stellar wind {v^ > Vorb), 
the Bondi-Hoyle accretion rate is 

Ma OC M^—^ 

(here a is the binary's semi-major axis). Substituting equation (|I]) for stel- 
lar wind mass loss rate and using M-L and M-R relations yield the X-ray 
luminosity 

~ Ma ~ M^'Va^ . 

Observations also suggest that the binary systems semi-major axis distribu- 
tion has a flat shape dXoga = const (Massevich, Tutukov 1988) and thus is 
independent of the optical star mass, so Mo ~ lI^^, dMo/dL^ ~ Lx^^^. 
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Then for the luminosity function of X-ray sources accreting from stellar wind, 



dN dN dM„ 



dLx dMo dLx 

using the Slapetdr mass function for main sequence stars dN/dMo ~ M~^ '^^, 
we get 

^ ~L-o-94L-3/5 ~L-i-5^ (4) 

IT a; X X \ J 

dLx 

which is very close to what is actually observed. In fact, in the stationary 
case, instead of using the Salpeter initial mass funciton, we should utilize the 
steady-state mass distrubution with account of the mass loss rate dMo/ dt ~ 
Muj, which would make the mass dependence even steeper than the initial 
mass function has, so the power-law index in equation changes towards 
higher values, as we need. 

In the case of slow stellar wind < Vorb, the rate of mass capture by the 
compact star is mainly determined by its orbital velocity, which for high- 
mass X-ray binaries with Mo/M^ » 1 depends on the optical star mass in 
the similar way as the escape velocity does, so no qualitative changes are 
expected. 

Note that although the chemical composition is important for absolute values 
of stellar luminositites and radii, it does not affect their functional depen- 
dence on mass which determines the shape of the luminosity function. 

Remarkably, a close to the obtained dependence can be formally derived in 
the case of accretion from optical star through the Roche lobe overflow, too. 
Indeed, in this case the accretion rate occurs in the thermal time scale of 
normal star, tkh- Mo ~ Mo/tkh ~ MoRqLo/GM^. Then again ~ Mo ~ 
Mo ~ 

_ r -0.94 r -0.64 ^ r-1.58 

J_y ryi -1—/ H-n J-V ^ 



dLx 



■'X -^x -^x 



However, in real high-mass X-ray binaries the mass transfer rate in the ther- 
mal time scale leads to the supercritical accretion regime, the accretion lu- 
minosity becomes on the order of the Eddington one and does not depend on 
the optical star mass. In reality, in stationary case it can be notably less than 
the Eddington value due to absorption in the powerful wind outflow from the 
accretion disc, as in the case of SS 433. Such sources are not numerous and 
can not significantly affect the shape of the X-ray luminosity function. 
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